function [U,V]=Lmafit_Demo(M_E,Omega0,r,x_size,y_size)

m  = x_size; n =y_size;

Known = Omega0;
data = M_E;
        
opts.tol = 1e-10;
opts.maxit = 3000;
opts.outobj = 1;

[U,V,~] = lmafit_mc_adp(m,n,r,Known,data,opts);
end